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Abstract 



A K-user pseudo-orthogonal (PO) randomly spread CDMA system, equiv- 
CN , alent to transmission over a subset of K' < K single-user Gaussian channels, 

is introduced. The high signal-to-noise ratio performance of the PO-CDMA is 
analyzed by rigorously deriving its asymptotic multiuser efficiency (AME) in 
the large system limit. Interestingly, the X'-optimized PO-CDMA transceiver 
Q I scheme yields an AME which is practically equal to 1 for system loads smaller 

than 0.1 and lower bounded by 1/4 for increasing loads. As opposed to the van- 
ishing efficiency of linear multiuser detectors, the derived efficiency is comparable 
^ I to the ultimate CDMA efficiency achieved for the intractable optimal multiuser 

«Q ' detector. 

o 
^. 

^ ■ 1 Introduction 

o 

Non-orthogonality of realistic spreading sequences, serving as signature codes in direct- 
Q ! sequence code division multiple-access (DS-CDMA), has an inherent detrimental effect 

in cellular communications. Nevertheless, a de-facto orthogonality between transmis- 
sions can be achieved via a proper joint binary signaling, for which multiple-access 
/\ • interference (MAI) agrees with information polarity. 

By borrowing analysis tools from statistical mechanics, counting the number of 
metastable states of the Hopfield model of neural networks [1 , we rigorously derive 
the asymptotic multiuser efficiency (AME, [2 ) of this pseudo-orthogonal modulation 
scheme for the common random spreading CDMA channel and draw its optimum value. 

2 Pseudo-Orthogonal Random CDMA 

Consider a perfectly power-controlled, synchronous, i^-users, A^-chips random binary 
spreading, binary signaling, Gaussian CDMA channel. A random CDMA channel is 

termed 'Pseudo- Orthogonal' (PO) if the transmitted signaling vectors { — 1, 1}^ 9 b = 
{bi, . . . jbx} are chosen such that 

Rb = Db, (1) 

or in scalar form 

K 

Y^p^^bi = Xkbk yk = l...K, (2) 



i=l 



where the K x K symmetric matrix R, with entries pki, is the spreading sequences' 
normahzed crosscorrelation matrix, and the K x K matrix 



^5 



D = diagjAi > 0, . . . , Xk' > 0, \k'+i > — oo, . . . , Xk > — oo} (3) 

is a diagonal matrix, in which the first K' < K diagonal entries must be positive and 
the rest can get any value. 

To put differently, the allowed signaling vectors b are those for which the MAI does 
not flip the information bits of a given fraction 7 — K'/K of the users. Evidently, this 
pseudo-orthogonality is gained at the cost of reducing the signaling entropy per user, 
H, to be less than unity (as opposed to multiuser precoding 3j in which there is no 
restriction on the input signaling.) 

On the other hand, this scheme reverts the multiuser channel to an equivalent 
set of 'jH single-user additive white Gaussian noise (AWGN) channels, with much 
simpler detectors w.r.t. the intractable optimal multiuser detector (MUD). This joint 
signaling scheme is particularly attractive for the CDMA downlink (but may also be 
utilized for the uplink in case of cooperation between users.) In order to evaluate the 
high signal-to-noise ratio (SNR) bit-error performance of the PO-CDMA scheme, its 
AME is analyzed. 

3 Asymptotic Multiuser Efficiency 

Theorem 1 In the large-system limit analysis, i.e., K,N -^ cxo, yet the system load 
factor (3 = K/N is kept constant, the PO-CDMA channel's AME, which characterizes 
the performance loss (in effective SNR) as the background noise vanishes '^21, is proved 
(for a given active users fraction •y) to get the form 

+ 7'log2(2Q(t*))+7(l-7), 

where Q{-) is the complementary cumulative distribution function of a standard Gaus- 
sian random variable^, and the fixed-points a* , b* and t* are devised by solving numer- 
ically the following set of equations 

/3-(l-i^)-. 7 W) . 



^Q{t) 



with an auxiliary variable t — {b — l)/y/a^. 

4 Proof 

If a particular multiple-access system achieves bit-error-rate (BER) per user P{(t) in 
the presence of MAI and AWGN with power spectral level equal to a"^, then the AME 
w.r.t. a single-user system (or a fully-orthogonal CDMA) is defined by j21 

r] = lim e((j). (4) 



lQ(x)^l/^/2^/;°dyexp(V/2) 



The term e{a) is the energy per user required to achieve BER equal to -P(o") in a 
single-user Gaussian channel with the same background noise level. 

According to definition (JH), the AME for the PO-CDMA case (P), ■r]{/3, 7), is nothing 
but jHl/S, 7), the fraction of 'orthogonal' users multiplied by their signaling entropy (as 
one is interested in an equivalent set of single-user channels with unit entropy.) Thus 
we have to compute the non-trivial entropy H{i3, 7) under the PO-CDMA transmission 
constraints (J2I). 

A i^-length binary transmission codeword b'^ = {b1, . . . ,fe^}, composed of all K 
users' bits at a given channel use, for which the PO constraints ^ hold, satisfies the 
condition 

/oo J^ /-oo ^ 

JJ dXkSC^ pkih - Xkhk) = / U d\k5(^ pkibi - Xkbk) = 1, (5) 

- fc=l i "^^1 fc=l i^k 

where the function 6{-) is the Dirac delta function. Condition ^ can be reformulated 



as 



/oo ^^ 
Y[dhSC^apkibi - aXkbk) = 
-1 fc=i i^fc 

/oo ^ 
Y[ dXk6(^ apkih - Xkbk) = 1, (6) 

'" fc=i ij^k 

where a = 1//3. 

Let the random variable N{f3, 7, R) denote the number of PO-CDMA allowable 
codewords, i.e. 

/oo ^ 
l[dXkS{J2c^Pk^b^-Xkbk). (7) 

-" fc=l 

r information rates, th( 
is defined ^, in bit information units, as 



fc=i j^fc 

Assuming equal user information rates, the corresponding asymptotic signaling entropy 



^ log2N(A7,R) .„. 

^(^7) = lim . (8) 

Assuming self- averaging property j3], in the large-system limit K ^ 00 the number of 
successful codewords N(/5, 7, R) is equal to its expectation w.r.t. the distribution of R, 
i.e. 



hm N(/?,7,R)=Ar(/?,7) 

1— >oo 

}"^ / \\'d^kY^\\\^\Z^^Pki\-^kbk^^ , 

'*°°-^-°'k=i b \ fc=i iT^k / R- 



where J\f{P,'y) and < ■ >r denote the average and averaging operation w.r.t. R, 
respectively, and J2h corresponds to a sum over all the 2^ possible values of b. 

Representing the delta function by the inverse Fourier transform of an exponent, 
expression (jH)) can be rewritten as 



/oo " -1 I'OO " 

-a r,_i V y J —oo r,_i 



fc=l ^ ' >'-""fc=l 



K' 



X ^ exp ij ^ uJkXkh) (exp(-j^ apkAuJk) ) , (10) 

b k=l ^ i^k / R- 

where j = -y/— 1 and u is the angular frequency of the Fourier transform. Substituting 

bkUJk for Uk, we find 



/OO 1 /"OO 

X ^ exp (^j J2 ^fc^fc) ■ E, (11) 



where 



E = / exp (^ - j J^ apkibibkUJk^ 
i+k 

exp ( - j ^ — ^ <sf feifofeO^fc j \ . (12) 



i-^k fJ.=l 

A 1 /ArV^A^ 



In the last equality the cross-correlations pki = l/^X]u=i ^k^i ^^^ expressed exphcitly 
as a function of the binary chips s^ = ±1. The expectation E can be also written as 



E = exp{ja^uJk)(exp(^-—J2(J2^khuJk){^s'^bk)j) . (13) 

k ^ fi k k ' ^ 

Using the transformation 

exp i-,AMK, ^ £ p^ £ ^^ ,14) 



/ K 7 7 

X exp (^jy (aj - 6j) - -^^^(a^ + b^) - -j=Bf,{a^ - b^ 



expression (fT^ becomes (here, and hereafter, logarithms are taken to base e) 

'•°° da., r^ T-r rf6„ 



X exp(^jy^(a2 -62) + ^log(cos(cfc,^))j, ^^5^ 

where 

Cfc,;. = ^(^fc(a^ + &^) + (a/,-6^)). (16) 

Since ^^ s^fe^ in (fT^ is 0(VK ) for a vast majority of codewords, for the expectation 

E to be finite, a^ and 6^ must be (9(1/Va). Hence, expanding the log(cos(-)) term 
in exponent ()15p and neglecting terms of order \jK and higher, we get 

e^pO"E-)/:n<^/:n "" 



E 



k 

K 
X exp I J 



(27r/i^)V2y_ii(27r/K)i/2 



k,fi 



where 

Ck^, = (ujlia^ + h^f + 2ujk{al - hi) + {a, - h^f). (18) 

Now, the solution of the iiT-dimensional integral |T7jl of the expectation E is per- 
formed using the following mathematical recipe: New variables are introduced 

« - ^E(«/^ + M', (19) 



M 



^ - ^Y.«-K) + ^- (20) 

Equations (fT^ and (j^U]) can be reformulated via the integral representation of a delta 
function using the corresponding angular frequencies A and B, respectively, 



oo 



dadA ,.^^,, v^('2u + ^u)^- 



exp (jirA(aa - }_^ ^^" ' ^' )) = 1, (21) 

Substituting these (unity) integrals into the expectation expression (fT7|) and rewrit- 
ing it using a and 6, the integrations over a^ and h^ are decoupled and can be performed 
easily. Next, for the asymptotics K —>■ oo, the integration over the frequencies A and 
B can be performed algebraically by the saddle-point method 13]. 

According to this method, the main contribution to the integral comes from values 
of A and B in the vicinity of the maximum of the exponent's argument. Finally, the 
E term boils down to 

/"°° da db 1 ^-^ 2 ^-^ 

1 {l-hY 1 
X exp(ira(6-- + ^^-^ + -loga)). (23) 

Substituting the expectation term (j^ back in (fTTj) . the integrand in the latter 
becomes independent of b, and therefore the Ylh '^^^ ^^ substituted by multiplying 
with the scalar 2^. Hence, 

^K-K' 



JJ d\k — ^^ / JJ duk exp \3 ^ t^fcA^ 
/•°° dadh /^ ,^ 1 (1-6)2 1, ,\ 

X exp (^- -aa^c^^+ja6^u;fcj, (24) 



fc fc 



where the resulting lo dependent integrand is a Gaussian function. Thus performing 
Gaussian integration and exploiting the symmetry in the i^-dimensional space, we get 



Kr<o ^ V 2^""^' r dadh / ^ ,^ 1 (1-6)2 1 

X expfK'logf./5rrfAexp(-M±^))). (25) 

V Wj aa J_Q, 2aa / / 



Using the rescaling {ab + \)/y/aa -^ A, the integral (j^ becomes 



AfiP, 7) = hm r /^ exp {Kg{a, b, 7, /?)), (26) 

where the function g[a, b, 7, P) is defined by 

^(a,&,7,/5) = ^(^'1+ ^^^ + ^loga) +7log(2Q(t)) + (1 -7)log(2), (27) 



with an auxihary variable 

A 



t ^ ^^^-^\ (28) 



Again, ioi K ^ 00, the double integral in (j26j) can be evaluated by the saddle-point 
method. Hence, we find^ 

U{(3, 7) oc lim exp {Kg{a\ b\ 7, /5)), (29) 

where a* and b* are found by the saddle-point conditions, which yield the following 
equations 

?il^ . r'(<i^-l)+7<^^0. ,30) 

da a Q[t) 



db "^ ' a ' Va^Q(t) 

The operator Q' denotes a derivative of Q w.r.t. its argument. This set of saddle-point 
equations can be solved numerically to obtain its fixed-points a*, b* and t*. 

Finally, substituting (j^ into (jHI the asymptotic entropy, in nats, is now easily 
obtained 

H{P,^) = g{a*,b*,j,P) 

= ^(b*-l + ^-^^p^ + 1 iog«*) + 7 log imn) + (1 - 7) iog(2),(32) 

which, along with equations (l ^ - ffHUjI -fl ^ . concludes our proof and forms the desired 
theorem. 



5 Results 

Fig. [T] presents the optimum AME of PO-CDMA ?7opt(/^) ~ ^'^^'^■yViPyl)^ drawn via 
an exhaustive search over all possible values of the active users fraction 7. Also drawn 
for comparison are the AMEs of the single-user matched filter (SUMF), decorrelator, 
and linear minimum mean square error (LMMSE) detectors [2^. 

Interestingly, ?7opt(/5 ^ 0-1) is practically equal to 1, which is the optimum AME 

of randomly spread CDMA [Zj, obtained for the optimal MUD. For large f3 loads the 
optimum AME of the examined scheme converges to 1/4, and not as for the other 
sub-optimal detectors. The corresponding optimal fraction 7 is found to decrease from 
1 to 1/2 as we increase /3. 

^Thc exponent pre- factor in (|29|) is not required for computing the asymptotic entropy, and there- 
fore it is omitted. 
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Figure 1: Optimum AME of PO-CDMA compared to optimal MUD and linear detec- 
tors. 
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